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A GENERAL KERNEL-BASED REGULARIZATION METHOD FOR COUPLED ELLIPTIC 
SINE-CORDON EQUATIONS WITH GEVREY REGULARITY 

VO ANH KHOA*, MAI THANH NHAT TRUONG, AND NGUYEN HO MINH DUY 


Abstract. Developments in numerical methods for problems governed by nonlinear partial differential equations under¬ 
pin simulations in diverse areas. This paper is devoted to the regularization method for the coupled elliptic sine-Gordon 
equations. The system of equations originates from the static case of the coupled hyperbolic sine-Gordon equations 
modeling the coupled Josephson junctions in superconductivity, and so far it addresses the Josephson 7r-junctions. In 
general, it is ill-posed due to instability of solution. Using the modified method, we obtain stability and convergence 
results. The results are observed as the generalization of many previous works. 


1. Introduction 


The Josephson junction is a quantum mechanical device, which is formed by two superconducting electrodes separated 
by a very thin insulating barrier. The Josephson 7r-junction is a specific example of that when no external current or 
magnetic field is applied. It is the corner junction made of yttrium barium copper oxide (d-wave high-temperature 
superconductivity). The fundamental equations modeling such a junction (see Chen et al. read 


du 2 EV du 
dt P ' dx 



du 

dy 



H 


1) 


J {x,y) = -Jo {x,y)smu, 

where u is functional to describe the relative phase between the superconducting metals I and /J, causing the 
Josephson tunneling current J per unit area which depends on the properties Jq of the barrier. Experimentally, 
E corresponds to the electron charge, V is considered as a time-and-space-dependent potential difference across the 
barrier, P is a Planck’s constant, d is a constant with respect to the London penetration depths for the metals, oj 
performs the speed of light, and Hj represents the x and y component of the magnetic field. 

After substituting those equations into the Maxwell equation, it yields the barrier equation or the hyperbolic sine- 
Gordon equation 


(l.I) 


d'^u 


d^u d'^u 


= —7 ^ sinu. 


where cq depends only on uj and 7 includes lo, P, E^d and Jq. 

The static version of dm thus reads 

d^u d'^u 
dx"^ dy'^ 

and if we proceed the change of variable v {x,y) = tt — u {x, y) under the homogeneous Neumann boundary condition, 
we are concerned with the equation 

d'^u d'^u _2 
-^ + = 1 smu, 

dx'^ dy^ 

which addresses us to the relation of the Josephson 7r-function. Mathematically speaking, the Neumann condition 
guarantees the equivalence between those two kinds of junctions, but the others, for example, the homogeneous Dirichlet 
boundary condition. 

From that motivation, the purpose of this paper is to consider, as the dynamical model for the coupled Josephson 
junctions, the following system in two-dimensional: 


= —7 ^ sinu. 


( 1 . 2 ) 


oil'll d'^u 

+ ai J- 71 sin (JiiM J- 5i2v) + anuP ai2V 


h, 
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(1.3) 


+72 sin ((521M + (5221’) + 0 - 2 lU + a 22 V = f 2 , 


where Ui € R+, 7 ^, , ( 7 ^ G K are physical constants and fi are called forcing functions, i,j £ { 1 , 2 }. 

There are hundreds of studies on the hyperbolic version of For instance, Levi et al. nni considered the 

chaotic behaviors of numerical solutions for the coupled hyperbolic sine-Gordon equations with damped terms. The 
more general system can be regarded as the problem of identification the physical constants by Ha and Nakagiri in 
[ 6 ] and the answer for the question of necessary condition to gain optimality for those parameters can be found there. 
Up-to-date, the amount of papers which are related to those problems still increases without cease. We, however, stress 
that the results as far as we know for the elliptic sine-Gordon equations are very scarce. In particular, we only hnd some 
theoretical and numerical researches, e.g. in [ 10 IIU 1 I 3 ] and references therein. Additionally, it should be mentioned 
that the Gauchy problem for such equations is ill-posed in principle in the sense of Hadamard where the stability of 
solution fails. Consequently, our analysis will underscore the so-called regularization method. 

In recent years, various types of regularization methods have been of continuous interest to researchers in a wide 
range of disciplines. We refer the reader to many interesting results in HI H n [3 n US- We notably remark the 
studies of Tuan and his co-authors on the nonlinear elliptic equations for which they are strongly related to our tackled 
equation (ll.2l) - (ll.3l) . In fact, the authors considered in [1^ the modified method to show that there exists a filtering 
kernel stabilizing the exponential instability when the spectral representation of solution is used. Only a short time 
ago, taking a more general nonlinear source term into account, they successfully obtained in [15] a new landmark by 
using the truncation approach. In light of the aforementioned works, we put ourselves into the study of the modified 
method to regularize the system (HID-(ESI)- 

Let (0, a) X (0,6) for a, 6 > 0, a couple of real unknown functions {u,v) is sought for {x,y) £ [0,a] x [0,6]. The 
problem given by (HISD-iHISD naturally along with the zero Neumann conditions: 


(1.4) 


du , . du , dv , , dv , 


and associated with initial conditions 


(1.5) 


{u (0, y ), V (0, y)) = (mq [v) , -^0 (y)), 


(1^ (0, y) > y)) = (wi (y) - ^1 (y)) > 


y e (0,6). 


In this paper, we consider problem 11.2111.51 and prove stability and convergence results of approximate solutions 
constructed by the generalization of modified method. 

2. A THEORETICAL FRAMEWORK OF REGULARIZATION METHODS BY SPECTRAL THEORY 

2.1. Abstract settings. Let (•, •) be either the scalar product in or the dual pairing of a continuous linear functional 
and an element of a functional space. The notation jl-jl stands for the norm in and we denote by INI;, the norm in 
the Banach space X. We call X' the dual space of X and denote by G ([0, a] ;X) for the Banach space of real functions 
M : (0,a) —>■ A measurable, such that 

ll«llc([0.a];A) = SUp ||m (x, •) || x < OO. 

0<a:<a 

We define the usual inner product and norm of (0, 6) by 

pb 

{u,v)= u{y)v{y)dy, ||m|| = (m, Vm, w £ (0, 6). 

Jq 

Let us also define the functional space 

Y = {u € (0,6) : Uy (0) = Uy (6) = O} , 

the closed subspace of (0, 6). 

It is significant to remark that £ Vfl C°° ([0, 6]) is an orthonormal basis generated by the operator in 

(0, 6) and it is associated with the eigenvalue which is such that 

0 < Al < A 2 < ... lim A„ = 00 . 


We turn to introduce the abstract Gevrey class of functions of order s > 0 and index 12 > 0, denoted by G®, defined 


by 


\ u€L^{0,b):^ |(n,(^„)|' < (xd , 


2 



equipped with the norm 


1/2 


\{uAn)f j < OO. 

We now make the following assumptions: 

(Ai) uq, Vo, ui,vi € (0, b) and the functions Ug, Vq, u\,v\ S (0, h) are measurement datum with noise level e > 0 

such that 

||Mo-Woll<e, lko-Go||<e, ||Mi-Mi||<e, Hiii - uj|| < e. 

(A 2 ) /i,/2eC([0,a];L2 

Our objective in this paper is to develop the computational foundations, we further assume that the system (1121)- 
(jl.5l) has a unique solution in O. In the next subsection, we shall give the regularization method and show our main 
theoretical results. 


2.2. Regularization method, well-posedness and convergence. Let {u,v) be a pair of solutions to the problem 
(1121)-dLll). So, the Fourier series corresponding to u and v are given by 

00 00 

u (x) = ^ {u [x] , (j)n) 4>n, v{x) = '^ {v (x) , (/)„) (/)„. 

n—1 n—1 

Multiplying each equation of the problem by (pn (y), and integrating with respect to y from 0 to b, we, after some 
computations, reduce the following system of Cauchy problem: 


(2.1) 

d^ 

{u {x) , 4>n) - aiXn {u (x) , 4)n) = {J^l (x, U, v) , (j)n) , 

(2.2) 

{U (0),())„) = (Uo,0n) , 

{Ux (0) ,())„) = {Ui,4>n) , 

(2.3) 

{V (x) , (j)n) - a 2 X„ {V (x) , (j)n) = {J^2 (x, U, v) , (pn) , 

(2.4) 

{V (0),())„) = (Ug, ())„), 

{Vx (0) ,(/i„) = («!>„) , 


where the functions J^i and J ^2 are given by 
(2.5) J'l {u, v) = fi- 7 i sin {SnU + Si 2 v) - cthM - (J 12 V, 


( 2 . 6 ) 


J ^2 {u, v) = f 2 — 72 sin {S 21 U + S 22 V) — (J 21 U — (T 22 G. 


By solving the problems (I2.1I) - (I2.2D and (I2.3I) - (I2.4I) . respectively, we say, under the nonlinear spectral theory, that u 
and V belonging to C ([0, a ]; (0, &)) are mild solutions to the nonlinear problem (ll.2l) - (ll.5l) if they satisfy the integral 

equations 


(2.7) 




n=l 


( 2 . 8 ) 


' (a^) = 


n=l 


cosh aiA„a:^ {uo, (j)n) + 
cosh (^^/a^x'j {vo, 4>n) + 


sinh (^y/aiXnx) 
\/ CHiXn 

sinh [\/a 2 Xnx) 
\/ OL2Xn 


(ui,(/)„) -f / 

Jo 

{Vi,(j)n) + 

JO 


sinh (\/aiA„ (x - ^)) 
\/ OilXn 

sinh (Va 2 X„ {x - ^)) 
V Q:2A„ 


{Tl (^, U, v) , (j)n) 


{T 2 {^,U,v),(j)n)d^ 


Therefore, let us observe (12.71) and (12.81) that when n —>■ 00 , the rapid escalation of cosh(a;) and tells us the 

ill-posedness of the considered problem in L^. A small perturbation in the data, in general, described by the assumption 
(Ai) can arbitrarily show a huge error in the solution. Thus, performing any computation is impossible in this case 
and then a regularization method is required. 

Now given /3 = /3 (e), let us define the filtering function (/3, x) : (0,1) x [0, a] —K for n G N, /c > 1, i G {1, 2} 

which is such that 


(2.9) 


(/3,a;) 


g-VaiA„(a-x) 

2/3-k ' 


Prior to investigating the use of our filtering function, we prove its boundedness by the following lemmas. 
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Lemma 1. Given /3 > 0, fc > 1 satisfy > kjd then for x G [0, a] the following inequality holds 


( 2 . 10 ) 


(^,x) <-^ (fca) “ /3 M ’ Vn G N,i G {1,2}. 


Proof. Let us first put g {z) = , then taking its derivative implies 

, , , Re~^^ — I3kz^~^ 

9 \^) =-9“- 

{/3z* + e~^^) 

We realize that g attains its maximum at z = zq where zq solves the equation e~^^° = ^Zq~^. It follows that 

1 


( 2 . 11 ) 


giz) < 


I3zt 


k0^k—1' 


0 i? ^0 


Moreover, by using the fact that > Rzq we have 


which leads to 


_ ^k-l Rzo < _J_ 


1, / i? 
•2^0 > TT; in 77T 


,kRzo 


kR 


kp 


Combining this with (12.111) . we conclude that 


, X 1 1 

9{z)<—^< 


kR 


Pzl P Vln(^) 

As a result, we go further to obtain the fact that for 0 < r < ii 


( 2 . 12 ) 


e-^^g{z) = 


{Pz^ 


o—Rz\ R 


{Pz^ 


o—Rz 


< {g{x)f ^ < {kR) 


r) 


p^ 


Rk\\ 


Therefore, the proof of the lemma is now straightforward due to replacing z = y/atXn, r = a — x and i? = a in 

dnn). □ 

Remark 2. In lieu of considering r = a — x in the last step of the proof above, we further take r = a — x + ^ for x > ^. 
Thus, we obtain 


(2.13) 


^n.fc,z {P,X- C) = 


^-\/ai\n{a-x+G) I 


k(.x-0 ^ 


2/3v'af A* + 2 


< t: (fca) “ “ In — 




kp J 


fG{l,2}. 


So now, we define a pair of regularized solutions under construction of the filtering function in such a way that the 
solutions do uniquely exist and obtain themselves the stability and convergence results. Here they are given by 


w = E 

C 

(2.14) 

n- 

oo 

v^{x) = Y. 


n—1 
oo 


n=l 


4 'n.fc,l {P,x) + 

1 


c — y/ai\riX 


(ug, pn) + 


1 


V ctiX„ 


'f'n.fc.l {P,x) - 


^-y/ai\nX 


{ui, pn) 




n—1 


4'„,fc.2 {P,x) + 


'I'n.fe,! {P,X - C) - 
~V 


c-Votl>^n{x-0 


{Vo^Pn) + 


\/ 0:2 Xn 


{Rl (^,U^ V*") , pn) d^pn 
g—\/Q2^ 

4'n./c,2 {P,X) --- 


{vl,Pn) 


(2.15) 


+ 


00 px -| 

y / 

Va2Xn 


4'n.fe.2 (/3,X - ^) - 


g-%/a7W(a:-{) 


{T2{^,u\v^) ,pn) df,pn. 


We notably mention that if we fix G C ([0, a]; Lp (0, 6)) in (12.141) . then we obtain the integral equation u'^ (x) = 
Q (u*^) (x) where Q dehned by the right-hand side of (12.141) maps from C ([0, a]; Lp (0, 6)) into itself. The boundedness 
of the hltering function by (12.101) and (I2.13|) in combination with the global Lipschitz of the functional yields the 
existence and uniqueness of the solution for every /3 (e) > 0 by the well-known Banach fixed-point theorem. The details 
can be directly followed in m- In the same vein, we also obtain the existence and uniqueness of solution to the equation 

dnsD. 
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Let us now consider the stability of (12.141) - (12.151) . 


Theorem 3. If and {U^, V^) are two pairs of solutions to \2.14\ l- li2.15\} . respectively, corresponding to the initial 

states {uq,uI,Vq,vI) and (UQ,Ul,VQ,Vf), then for all 0 < x < a, 


K(x)-C/^(x)f+ ||z;^(a;)-y^(a;)f “ E {\\< - + U - Vff) , 

where Ca is a positive constant only depending on a,ai,a 2 , 311 , 612 , 621 , 622 , o'ii,ai 2 ,cr 2 i, cr 22 - 

Proof. From (I2.14I) - (I2.15|) . one deduces that 


\{U^ (x) - {x) ,(j)n)\ < \'^n,k,l{l3,x) + 


^—y/aiXnX 


yj G^i A 

. 

7o \/ CklAr: 


d'n.fc,! (/3,X) - 




'^n,k.l {P,X- f) - 


|(uo {7 q,0„)| 

\{u\-Ul,cf^)\ 

I (j-i (e, C/% </>„)! de, 


3-Vo!iA„(a:-{) 


\{V^ {x) -V^ {x) ,(j)n)\ < {'^n,k, 2 iP,x) + 


a —\/a2 A„x 


+ 


1 


y/af^Xr 

X’ 


^'n,fe ,2 (/3,a;) - 


OL-2XnX 


1 


^n.fc,2 {I 6 ,X- f) - 


- yQXn)\ 

|(uj - Vf,(j)n)\ 

\{P2Xu^,v^)-XXU\V^),<fn)\df. 


^-Vap\Z{x-^) 


10 \/ 02^. 

Since we have, by (I2.5F(I2.6F that 

(2.16) \Ti (e,u%u% - Ti X C/%y%| < (I 71 I Idiil + kill) - C/% + (I 71 I ki2| + ki2|) , 


(2.17) 


\T2 Xu^,V^) - T 2 X UWX < (I 72 I k 2 l| + k 2 l|) - C/% + (I 72 I Ifel + k 22 |) X - ^^1 , 


they thus give us, by Parseval’s relation in accordance with the boundedness of the filtering function and the very 
basic inequality (a + 6)^ < 2 {of + 6^), that 



+ 1 1 ( Ikn - C/n%l V Ikg - Pol'" 


+ 1) + ^-^1%!") 


(2.18) 


- ({ka) 

Jq ckAi \ 



2kU-x) 


+ lj diOdf, 


where (i(a;) = k)ll + Ik*" (2^) - (a^)ll i« = min {ai, 02} > 0, C = ^£{1,2} (17% Ikjl + ky |) . 

> 1 for all X G [0,a], it then follows 


2kx 2x ( ( ^ \ "N 

Assume that /3 S (0,1) is small enough such that (ka) “ /3“~ (In ( 1 1 

from (12.151) that 



ko - + Iko - k.%1' + ^ (ik) - Xf + Ilk - Vff) 


2kU-x) 


(2.19) 


d{Odf. 




















































2fcar 2kx 

Multiplying both sides of (12.191) by {ka)~~ (3^ ^In ” and putting w (x) = {ka)~~ /3^ ^In " d (x) 

and thanks to Gronwall’s inequality, we obtain 


w {x) < 


ll«o - U^f + Iko - ^ (hi - U!f + - v,^f) 


2Cx\ 


or it can be written as 

( 2 . 20 ) 

_ 2kx 

2kx 2x ( ( C^\ 

d{x)<{ka)^r^ 


ho - c^of + ho - ^0 f + ^ (h) - uif + ||uj - y/f) 


exp 


f2Cx 


y CtAi 


This completes the proof of the theorem. 


□ 


In order to prove the convergence, we follow the strategy that proves the convergence between the exact solutions u 
given by and the regularized solutions corresponding the exact datum. Let us consider the following theorem. 


Theorem 4. Let us define {U'^,V^) by 




( 2 . 21 ) 


n=l 




^ — y/OLlXnX 


(uq, 4>n) 


'J ciiA„ 




3 ^iXjiX 


{ui , 0n) 


oo ^ 

s/ 

_1 ^0 


^ Jo 'J ciiA„ 


{P, X-S,)- - --- (J'l (C, V^) , 4>n) d^fin 


vhh = E 


( 2 . 22 ) 


'^n,k ,2 (/3, x) + 


o — \/Oi 2 XnX 


('^0: 0n) 


V Q!2A„ 


^'n.fe.2 (/3,X) - 


O—t/c(2A„X 




OO ^ 


^ 7o 


1 / g-Va2A„(x-^) \ 

^n.fc.2 (/?,X-0--- ,Pn) difin 


Suppose that the solutions given by h2.l\) - h2.8[) are such that u G C ([0, a]; G^)) , G C ([0, a]; GJ^) and v G 
C ([0,a] ;G®^) ,Vx € C ([0,a] ;G®^) for i^i > ^^,1^2 > and si = k, S 2 = k — 1. Then for fi G (0,1) the error 
estimate over (0,6) between {u,v) and {fA’^,V'^) is uniformly given by 

_ 2kx 

||u(x) -W{x)\\^ + ||u(a;) - V"(a:)f < Ca{ka)^ ^^(i-f) 

where Ca is a positive constant only depending on a, k, Oi, 02, Ai, 71,72, (5ii, 612,621, <522: cin, cri2, o’2i, 0 ' 22 - 
Proof. Observe it is easy to verify that 

(u± , (jjri} 


/ 


-Vo!iA„{ 


(u(x),<^„) + ^^h(^)^ = e^h(uo,<^„) . ' I 

x/aiAn \ VOllAn JQ 


{Pi {f, U, v) , fin) d^ , 




{v{x),fin) + YYYA=e'^^^^ ({V 0 ,fin)+ -^ ^ 

y OL^An \ yOt^An JO yOi^Ar, 


I 


^ p —\/a2An^ 


{P 2 if,, U, v) , fin) df . 


We put di (x) =u{x) — W {x) and d 2 (x) = v (x) — (x). First, we have the following 


di {x) = E 


P'JoJlXi 


2 PiJa\\i + 2 e-'^ 

1 

0 JoiA 


{U (x) , fin) + 


{Ux {x) , fin) 

'/oiXn 


■E 

n=l 


g-VaiA„(x-{) \ 

'^n,k,l{P,X- f) --- {Pi {^,U,V) - Pi {£,,W,V^) ,fin) dffin 




























































Therefore, by Parseval’s relation and (12.911 together with (12.101) . (12.131) and (12.161) . we get 


IMi(a;)||^ < 2/3^ ^ {u {x), 4)^) + \J a\ ^\n {u^{x) ,(t)n) 

n=l 


+2a^V / 


n JO 
n—1 ^ 


2kx 2x 


'^n,k,l {P,X- C) - 


„-VQiA„(a;-5) 




< (3 (ka) “ /3 “ (In ( ^ 


OO 


n=l 


(2.23) 





di (Of+ ||d2(C)ind? 


where Ci = (|7i| |(5ii| + \aii\f + (|7i| |(5i2| + |cri2|)^. 

In the same technicalities, ^2 (a:) can be estimated as follows: 


2kx 2x 


\\d 2 (a:)II < /3 (ka) “ /3 “ In 


(2.24) 


02^1 


“2 (^2 ||r;||^(-[(j^^].Q»i^ + ||a;a;||(^^[Q 

^ (fca) ‘ “ ’ /3^ (^In (lMi(Of+ IM2(0f)c?0 


where C2 = (I72I J Oil + 1 ^ 21 1 )^ + (I72I \S22\ + |cr22|)^. 
Combining (12.231) and (12.241) gives us the fact that 


w 


2 a^ (Cl C 2 


(x) < (a^-i + ^ pi + ^ / w{Od(. 

Ai \ai 02 ' ' 


where we have putted 


{x) = {ka) “ ^ (in (||di(Of+ ||d 2 (a:)f) : 


Ca — ai ||M|lc([0,a];G;() + ll“® II C([0.o];G“2 ) + “2 1^11 c( [0,a] iGji ) + f ®llc([0.o];Gj2) • 

By Gronwall’s inequality, we thus obtain 

_ 2kx 

,fe\\ a /2a^x (Cl C 2 


110 (a:)f + 11^2 (a:)f < Ca (oi ^+02 ^){ka) “ 0 (In 

which leads to the proof of the theorem. 


exp 


Al V Oil (X2 


H- + a — X 


□ 


Theorem 5. If P = e™ for m S (0,1], then the error estimate over Lp (0,6) between the exact solutions {u,v) given by 
and the regularized solutions {u‘^,v'^) given by ^2.14^ - i2.15\} is 

_ 2kx 

||u (x) - (Of + Ik (0 - (Of < Ca (kaf-^ (e2-(i-f) + 0(1-^)) (^In “ , 

for alio < X < a, where Ca is a positive constant only depending on a, fc, oi, 0:2, Ai, 71,72, dii, ( 5 i 2 , Oi, 622, cth, cri2, (T21, <722 

Proof. The proof is straightforward by Theorem [ 3 ] 13 ] and (Ai). In fact, choosing P = e™,TO G ( 0 , 1 ] we have 
Ik (re) - u'^ (:r)f + |k (x) - C (x)f < 2 (^\\u (x) - W (x)f + \\W (x) - u'^ (x)f + |k (0 “ V" (x)f + || V" (x) - v^ (x)f 


< Ca(fca)"^ 

This completes the proof of the theorem. 


a 

fee" 


□ 

































3. Concluding remarks 


We have successfully extended the modified method to the system of elliptic sine-Gordon equations wherein the 
Gevrey regularity plays a strongly powerful position in our analysis. The method is shown to be good convergence in 
comparison with well-known methods, such as quasi-reversibility method and quasi-boundary value method. In fact, 
one of superficial advantages is described by the error estimate at a; = a is valid in logarithmic type, i.e. 


||r (a) - (a)f + ||i; (a) - (a)f < C. (l + )) (^In < C. 

which cannot happen in the study of quasi-reversibility method (|9]). 

Moreover, the convergence rate is quite general. In principle, it is of the order O 
generalization of many previous results, e.g. oinj. 
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